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The  uniaxial  cornpfCssive  niech^^  propellant  is  well  documented 

over  the  temperature  range  from  -40  to  60  degrees  Celsius,  and  strain  rates  from  quasistatic  to  IR* 
sec-'  using  drop  weight  (Lieb  m^^  Hopkinson  bar  (Lieb  et  al  19119),  and  servohydtaulic 

(Gaxbnas  1991;  Gazonas  and  Ford  1992)  test  apparatuses;  Uniaxial  compressioa  of  rightcirculm 

cylinders  Of  M3 0  induces  fracture  damage  consisting  of  axial  cracks  that  grow  and  eventually  coa¬ 
lesce  to  form  macroscopic  conjugate  shear  fractures  at  large  strains.  MicrOgraphic  evidence  reveals 
that  microcracks  initially  form  between  the  subaxially  oriented  nitroguanidine  crystallites  (energetic 

filler)  and  the  nitrocellulose  (binder).  The  stress-strain  response  of  M30  is  slightly  nonlinear  prior  to 
fhe  rtlaxirtium  stress  level  (failure  stress);  the  propellantsubsequentiy  woriusoftens  until  ultiniate^^^^^^^^^^^^^^^^^ 

failure  or  rupture  occurs  (Gazonas  and  Ford  1992). 

^ 

studies  relatecl  to  their  constitutive  characterization  are  scarce.  A  recent  study  characterizes  the  ^ 
viscoelastic  response  Of  M30  propellant  and  shows  that  the  propellant  exhibits  nonlinear  (strain-^^^^^^^^^^^^^^  ^^^ 

dependent)  power-law  relaxation  Over  the  time  interval  lO’-  to  10^  milliseconds  (Gazonas  1991).  A 

complete  constitutive  description  of  the  solid  propellant  phase  is  Critical  for  accurate  description  of 

combustion  in  numerical  models  (Gough  1990)  Of  the  interior  ballistic  (IB)  process  because  the  rate 
Of  mass  generation  of  the  gaseOus  phase  during  combustion  is  proportional  to  the  amount  of  exposed 
propellant  surface  area  (see  Military  Explosives  1955).  Early  models  of  Combustion  assume  that 
time-dependent  surface  area  is  only  a  function  of  the  differential  changes  in  the  initial  propellant 

geometry  caused  by  denagration.  These  models  do  not  account  for  an  increase  in  surface  area  due 
to  deformation  and  fracture  of  the  prOpellaht.  Subsequent  IB  numerical models  incorporate 
eftectOf  enhanced  mass  generatiOrt  iate  due  to  fracture  by  using  .Surface  area  "multipliers"  (Keller 
arid  HOrSt  1989)  It  is  anticipated  that  the  constitutive  and  damage  characterization  of  single-grains 


of  propellaht  will  provide  insight  into  the  physics  govemirig  the  bulk  defbrmM^ 

lant  beds.  ' 

This  paper  employs  a  uniaxial  specialization  of  a  general  three-dimensional  constitutive 
theory  for  viscoelastic  materials  with  damage  (Schapery  1981).  Several  features  of  the  constitutive 
theory make  it  an  attractive  candidate  for  modeling  the  constitutive  behavior  of  M30  propellant; 
First,  the  theory  can  predict  the  observed  wotk-softening  behavior  in  M30  under  monotonically 
increasing  deformation.  Nficrocrackihg  materials  that  exhibit  work-softening  behavior  pose  Special 
problerns  for  cbristitutive  ntodelers  since  it  can  be  argued  that  ''shear"  fracture  planes  that  develop  in 
many  raaterials  in  compression  are  structural  or  geometric  features  that  corrupt  detection  of  the  true 
material  response  of  the  material  as  it  work-softerts.  Additional  problems  associated  With  a  loss  of 
hyperbolicity  of  the  Wave  equation  in  damaged,  Work-softening  materials  have  recently  been 
addressed  using  a  nonlocal  elasticity  approach  (Valanis  199 1 ).  Secondly,  the  constitutive  equations 
can  be  transformed  to  those  of  nonlinear  elastic  materials  through  correspbndence  principles;  The 
transformation  facilitates  the  solution  of  boundary  value  problems  encountered  in  the  theory  of 
nonlinear  Viscoelasticity .  Thirdly ,  micrOcracking  in  M30  is  characterized  with  a  damage  function 
that  is  related  tO  tirne-dependent  surface  area  evolution  in  the  pfopellant.  This  relation  could  be 
incorporated  into  IB  numerical  codes  for  more  accurate prediction  of  surface  area  evolution  and 
mass  generation  rate  during  propellant  combustion.  Finally,  the  theory  is  general  enough  to  success¬ 
fully  describe  the  nonlinear  viscoelastic  response  Of  a  variety  of  other  inaterials  that  include  :  marine 
sediment  (Schapery  and  Riggins  1982),  rocket  propellant  (Schapery  1982),  and  ice  (Harper  1986; 
Harper  1989).  Thus,  material  constitution  and  damage  evolution  in  a  variety  of  materials  can  be 
compared  within  the  framework  of  a  single  theory. 

Even  though  "material"  damage  is  treated  herein  as  a  scalar-valued  quantity,  predictions  of 
stress  versus  time,  failure  stress  versus  failure  time,  and  failure  stress  versus  strain  rhte  quantitatively 
agree  with  observations  from  isOthCrmalv  uniaxial,  cOnsttmt  strain  rate  compression  tests  on  the 
propellaiit.  It  isshown  that  microCracking  in  M30  is  characterized  by  an  “ellipsoidal”  damage^^^^  ^^^^^^^^^^^^^  ^^ 


furtctionr  iwivied  troni  shifted  isothermal  constant  slxairt^^r  data.  The  foritv  of  the  OlUpsoidal 

ilamage  function  is  compared  With  damage  functions  that  have  been  developed  for  0 

such  asv  ice  in  compression  (Harper  1986),  and  oil-shale  under  dynamic  blast  conditions  (Grady  and 

:-::-^:(:::WiZ;'THE:'tbNSTlTUTlVETHEORY::i;;:::iv'-:i: 


The  uniaxial  nbnlineaf  viscoelastic  cOhstitutive  equation  for  materials  that  possess  a  random 


or  regular  distribution  of  microcracks  can  be  Written  with  a  so-called  “modified  superposition  in^ 


graf  ’  (Schapery  198  l  ;  Schapery  1989).  The  uniaxial  strain  e  in  a  material  subjected  to  a  uniaxial 

stresS'Ohan  bevwTitten:as(^::y'V 


e(t);;;=;:Ef;..T;'ixt:-:t;)^  ,  ffffff:e;f.;':b;:(i) 


The  integral  in  (1)  is  also  know  hereditary  integral  or  as  a  convolution  of  the  functions  D  and 

f  Material  nonlinearities  and  damage  are  usually  incOrpo 

t(ct.S„)  .  (2) 


in  (1),  E^  is  ah  arbitrary  Constant  reterred  to  as  the  “reference”  modulus  with  units  of  stress.  Elastic 
behavior  is  ob  tained  When  D  ^  E^^.  D(t)  is  the  linear  viscoelastic  creep  compliahce  if  all  material 

nonlinearities  are  ihcOrpOrated  into  f.  The  creep  Compliance  is  definedas  the  Strain  response  nor¬ 
malized  to  the  unit  stress  input  (he.,  D(t)=e(t)/cy,,,  with  a(t):±CT^H(t)v  and  HfO  is  the  Heavis 
function  defined  as,  HfO  =  l  forf' >  0,  and  H(0  =  0  for  C  <  0).^ 

damage  parameters  that  influence  the  time-dependent  strain  in  f  l)v  The  <y- subscript  in  (2)  refers  to 
damage  parameters  developed  for  stress-history  inputSv  The  damage  parameters  are  used  tor  ; 
Strain-history  inputs.  In  (2),  f  is  written  in  product  form  with  a  power-law  stress  function  gj  and  a 


damaj^6  fuilction 


f  = :  gi(o)  e  sgn{<t) 


where  dj*  r  and  X  are  positive  cbnstam  The  stress-history  deperident  damage  parameter  is  derived 
(Schapery  1981)  by  integrating  the  itlation  between  crack-tip  Velocity  add  the  J-integral  and  is  given 


j  ^  f  ^  f 

a  J  <T,  1 


where  d,  and  q  are  positive  constants,  and  I  I  denotes  the  absolute  value  of  the  quantity.  A  damage 
parameter  similar  to  (4)  is  derived  by  Wnuk  and  Kriz  (1985)  by  integration  of  the  "KachanOv" 
equation  which  relates  the  rate  of  damage  growth  to  a  poWer-law  function  of  the  net-section  sdess. 
The  functions  g  j  and  g^^,  and  the  above  constants  may  be  differeht  for  characterizing  material  behav- 
ior  in  compfession  versus  tension  or  fOr  unloading  after  significant  plasticity.  The  genem 
the  theory  can  be  verified  if  the  above  constants  and  functions,  gj  and  gj-,  are  unique  for  a  variety  of 
stress-history  inputs  in  (1).  The  sigrtum  function  in  (3)  is  defined  as  sgn(d)  ^  1  for  compression  in 
this  Studyv  f,  is  a  "crack-tip  material  Coefficient"  which  can  depend  on  time  and  temperature  and 
material  aging  effects  (Schapefy  1981).  In  the  present  study  f,  =  1 .  The  damage  function  gj^(S^)  in 
(3)  reflects  material  damage  due  to  microcracks.  The  exponential  form  of  was  originally  pro¬ 
posed  by  SChaperyf  1981)  and  was  later  used  to  model  microstructural  damage  in  ice  subjected  to 
uniakiai  compression  (Harper  1986).  Later,  it  is  Shown  that  the  exponential  damage  function  which 
characteriz:es  damage  in  ice  does  not  satisfactorily  characterize  damage  in  M30  propellant  Instead, 
ah  “ellipsoidal"  damage  function  is  utilized  that  is  directly  determined  from  the  test  data.  From  (3) 
we  See  that  in  the  abserice  of  damage,  S^  =  0  and  g2^  =  1,  and  (1)  then  predicts  strain  in  a  nonlinear 
viscoelastic  material  without  damage.  Material  symmetry  changes  due  to  darnage  induced 
ropy  are  not  addressed  in  this  paper  and  damage  is  treated  as  a  scalar-ValUed  qUantityi  More  discus- 
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r 


Sion  of  this  topic  and  hi 
Kjrajcihovic  (1987)  and 


in 


Weitsman  (1 98 8). 


Equations  (1)  through  (4)  are 


the 


papers  of 


and  a  rriaterial’s  strain 


suitable  for  charaCteriiring  damage 
response  if  stress  is  a  controlled  input  for  the  test.  However,  if  strain  is  a  controlled  input  for  the 
test,  then  these  equations  must  be  inverted  in  order  to  predict  stress  as  a  function  Of  strain  history 


Where  g^^fSp  -  and 


where 


o(t)  = 


<7,  0” 


1/r 


g-(SJsgn(r) 


the  pseudo-straih  e"  isf elated  to  strain  history  by 


-  Ef  J  E(t - 1 )  rfe  dx 
0 


E(t)  is  the  relaxation  modulus. 


dx 


(5) 


(6) 


that  the  utility  of  pseudo-strain,  as 


Schapery  ( 1 982)  points  out 

a  strain  measure,  lies  in  that  fact  that  stress  versus  pseudo-strain  plots  are  single- valued  or  “elastic 
like”  for  Cyclically-strained  maferia^^^  How-ever,  the  e--Se”  transformation  in  (6)  does  not  produce  a 


tution 


The 


of  (5)  into  the  time-derivative  of  (4). 


or 


S, 


damage  parameter  S  .  is  obtained  by  substi- 


J 


Rearrangement  and  integration  of  the  result  leads  to 


(7) 


t 


Je  tj 


0 


The  lower  limit  in  (7)  is  zero,  which  corresponds  to 


dt 


a  no-damage  condition  at  time  t 


exponential  form  g,^(S^)  -  is  used  in  (7),  then  the  damage  function  becomes 


g,Tsj^(i+; 


,-I/q 


(8) 


0.  If  an 


(9) 


The  derivatioh 


of  a  more  general  damage  function  is  obtained  if  one  assumes 
given  by  (3).  Substitution  of  the  first  time-derivative  of  (4)  into  that  of  (3)  results  in  a  honlinear^^^^^^^^^^^^^^^ 
differential  equation  (Bernoulli  equation)*  solvable  for  stress,  which  can  be  linearized  and  integrated 


a  product  form  for  f  as 


using  a  power-law  function  of  the  pseudo-snain  as 


f=e“ 


an  integrating  factor  (Appendix),  The  identity. 


(10) 


is  used  in  the  derivation  and  is  obtained  if  E  D  ^  1  .  where  E  and  D  are  the  Carson  transformed 


relaxation  modulus  and  cn 

is  then  given  as 

;ep  compliance  respectively. 

A  generalized  form  of  the  damage  function 

q  ^q/r 
V  r  e” 


Ifgaa(S„)==  e^"in(l2),thendg,^^^^^ 


2.1. 


5^4 


(12) 


rate  e  input, 


e(t) 


L  This  section  provides  expres- 
.  For  a  constant  strain 


e  t  H(t) 


(13) 


Substitution  of  (13)  intt>  (6)  relaxation  modulus  bf  M3()  d 

(Gazonas  1991);  E(t)  =  Ej  t  "^  <  n  <  1),  and  Ej  =  E);  yields  the 


6  (t)  =  e  t  /(I  -  n) 


Bxperiffierital  results  for  isothennal  uniaxial  corhpression  of  M30  at  four  stfain  rates  appear 
in  Figure  1.  Each  stress-tirne  curve  is  a  composite  curve  formed  from  the  average  of  five  tests.  The 
compression  tests  are  performed  Utilizing  a  servohydraulie  test  appafatns  (MTS  810  ffi 
System)  that  is  described  in  mgre  detail  elsewhere  (GaZonas  1991);  The  maximum  piston  velocity  is 
on  the  order  oiM2  meters/sec  (39.4  feet/sec)  (and  this  limits  the  axial  compressive  strain  rate  to  500 
sec‘‘  in  25.4  mm  (1  inch)  long  specimens.  Constant  strain  rate  tests  are  performed  by  computer^^^^  ^^^^^^^^ 
control  Of  the  piston  velocity  via  feedback  from  an  extemally  raourtted  transducer,  linear-variable 
differential  transformer  (LVDT),  MTS  Model  244.11.  Force  is  measured  with  a  60  kN  (13.5  X  10^^^^  ^ 


lb),  qiiaitz-piezoelectric  ferce  (lansducet'v  Kistler  is  mounted  bh  the 

rnbvihg  piston.  Specimen  displacements  are  corrected  for  apparatus  distortion  which  has  a  ttiea^ 
stiffness  of  about  91.9  kN/mm  (52.4  X  lO^  lb/in).  Specimen  stiffness  ranges  from  4.5  to  14.0  kiST/m  m 
over  the  sfrain  rates  investigated.  Tests  are  conducted  at  a  room  temperature  of  22  ±  1  degrees 
Celsius./^:,;. 

Right-circular  cylinders  of  M30  gun  propellant  (Radford  lot  #128  B)  are  prepared  by  cutting 
specimens  from  six-inch,  solid,  stick  propellant  using  an  Isomet,  doubIe-bladed  diamond  saw.  The 
inert  lubricant,  molybdenum  disulfide,  MoSj,  is  sparingly  applied  to  the  specimen  ends  to  reduce  end 
friction  effects  and  test  variability  (GazonaS  and  Ford  1992).  The  chemical  Composition  and  homi- 
nal  speCimen  dimensions  of  M30  appear  in  Table 


first  determining  constants,  (Jj,  n,  r,  ^^,  and  q,  and  the  damage  function  g^^  in  (14)  through  (16). 

The  material  constants  are  determined  from  a  limited  Set  Of  experimental  data  (constant  strain  rate 
tests)  for  prediction  of  material  response  under  more  general  input  histories. 

The  damage  function  and  constants  are  obtained  by  plotting  log,o  gj^fS^)  in  (16)  versus  log, ^ 
in  (15)^  Experimental  cr(t)  data  are  used  in  (16).  Daniage  curves  for  a  set  of  arbitrarily  chOsCn 
constants,  Cj,  r,  n,  and  q,  appear  in  Figure  2.  However,  if  twO  curves,  taken  at  the  strain  rate  ^ 
tremes  (.01  seC''  and  42()  Sec  ‘)  are  shifted,  a  “master”  damage  function  is  formed  With  an  appropri¬ 
ate  choice  of  constants  (determined  by  trial-and-error,  see  Figure  3a)  which  are:  cr^  670  MFa  (97.2 
X  10*^  pSi),  r  1.2,  n  =  0.1,  q  =i  9,  and  e,^(  =  .01seC'‘.  damage  function  takes  an  “ellip¬ 

soidal”  form  (solid  line  in  Figure  3b)  in  logarithmic  coordinates  and  represents  the  “softening”  effect 
of  the  microcracks; 


8 


L5(  1  -  I -^)  > 

and  the  miGrocrack  “saturation”  level  is  at  log,,,  (S^^)  -  10. 

herein  for  M30  propeUant  is  compared  to  da 

for  ice  is  based 

linear  and  proportional  to -1/q  (see  (9))  whereas  experimental  data  for 
(Figure  3a):  iTie  damage  function  for  dynamic 

;:::^:■■■:^)(^0;0;''■■lion■of■oil-shale,;;:O0■)0;0.::;0;■^^ 

-  :: ::  (18)::::; 

:0:;0'::'::-::''::'>renceon;strain''rate'0  ::■:■ 'V 

;:;;;;(;:  ::):':0:.;  =::B  e ;t(l-  D)::’:;;:.''0:';::o0  (■;.::,( 

Damage  is  defined  here  as  a  scalar  quantity,  0  <  d  <  Iv  and  d  =  0  corf^^^  to  a  no-damage 
condition,  whereas  D  =  1  corresponds  to  complete  materialfailure.  In  (18)^  a  and  m  are  constants, 
and  in  (19),  e  is  lhe  sdaih  rate,  t  is  time,  and  B  is  the  indinsic  elastic  modulusv^^^^^^^^ 


OQ 


(Harper  1986;  Grady  and  Kipp  1980).  The  ultimate  utility  of  the  damage  function  lies  in  its  ability 
to  accurately  model  material  behavior,  regardless  of  the  method  Used  in  its  determinatibn.  However, 
models  developed  from  microstructural  considerations  are  iritrihsically  more  appealing  than  those 
developed  from  empirical  data  fits.  ^  ^  ^  ^  ^ 

1  Significance  of  Parameter  D:  Fracture  Surface  Area.  A  number  of  authors  have  attrib¬ 
uted  microstructural  significance  to  D  which  is  assumed  to  be  proportional  to:  1)  the  ratio  Of  ^ 
microcrack  area  to  total  area,  A7A,  (Schapery  19^^  Lemaitre  1985),  2)  ratio  of  defect  density  to  a 
"saturation"  defect  density,  d/d^,  (Rousselier  1981),  3)  the  ratio  of  the  radius  of  a  single  spherical 
microcrack  to  the  volume  Of  a  representative  unit  cell,  aVV,  (Budiansky  and  O'Connell  1976),  and 

4)  statistical  distributions  Of  the  ratio  of  the  number  of  broken  bonds  in  ''bundle''  iTiOdels  to  the  total 
number  Of  bonds,  n/N,  (Bolotin  1969),  to  name  a  few.  For  reasons  of  relative  simplicity,  ^  ^  ^  ^  ^  ^ 


microstructural  models  are  developed  for  cracks  which  grow  normal  to  an  applied  tensile  stress  field. 
For  materials  in  compression,  mixed  mode  cfaCk  growth  and  interaction  considerably  complicates 


stress  intensity  factor  (COstin  1987).  It  is  not  the  intent  of  this  study  to  develop  a  microstructural 


model  for  microcrack  growth  in  MSO  in  compression^  However,  D  coUld  be  expressed  in  tenns  of 
the  ratio  of  micrOcrack  area  to  the  "saturation"  microctack  area,  A/A,  (motivated  by  the  form  of  g^^ 
in  (17)),  Figure  4b  compares  D  vs  Log- firae  for  M30,  ice.  The  time-rate  Of  change 

of  D  increases  with  strain  rate,  i.,e.  ifbM  =  f(^^ 

Observations  indicate  that  the  amount  of  damage  (fracture  surface-area)  in  M30  propellant  is 
relativeiy  insensitive  to  strain  rate  from  10^2  sec'i  to  100  sec  '.  Over  this  strain  rate  range,  fracture 
surface-area  produc  tion  is  primarily  dependent  on  the  amount  of  axial  specimen  strain  ( i.e.  D  =; 

(e) )  and  secondarily  dependent  On  the  deformation  temperature  (Gazonas  et  al  1991).  However, 


thedreticai  predictions 

laagmcnt  size  decrcisas  strain  iite  incmases^^C^^^ 

in  oil-shale  decreasesfixim  0.1  tosOl  meters  as  in  die  moderate  loading  rate 

regimefrom  10  sec-^  to  too  sec-*;  the  damage  functionahd  material  const^^^^^ 
study  may  not  appropriately  model  the  constitutive  response  of  MSO  propell^ 
loading  conditions.  Additional  material  data  obtained  at  large,  dynamic  Spain  rates 

accurately  predict  materialbehavior  over  the  wide  spectrum  of  strains  and  Ibadmg  rates 

by  the  propellant  in  the  gun  during  firing.  The  hextsection  compares  data  and  c 
tions  of  time-dependenrstmsseS  in  M30  in  order  to  verity:  the  theory  and  ilte 


Time  dependent  vStfesses  for  constant  strain  rate  deformation  Of  M30  are  predicted  by  solving 
(17)  for  g,^(S^)  and  substituting  the  result  into  (16).  Stress  versus  time  predictions  are  plotted  in 

Figure  5  (solid  lines)  for  the  strain  rate  extremes  at  which  the  damage  lunction  is  defined.  Experi¬ 
mental  data  (symbolslare  also  plotted  for  comparison  with  the  theoretical  predictions.  Superposed 

ctifves  at  intermediate  strain  rates  (0.9  sec-l  and  89:sec^ 

theory  since  g5j(S-)  is  determined  froin  data  shifted  from  the  Spain  rate  extremes.  TTte  theoretical 

expression  in  (16)  provides  a  good  approximation  to  actual  stress-time  datapbtained  at  constant 

Spain  rate:  A  Cursory inspection  of  the  superimposed  SPess  versus  time  curves  revealathat  the^^^^^^^^^^^^^^^^ 

theory  accurately  predicts  maximum  Or ‘‘failure"  sPesses^  teilure  times,  and  the  work-softenmg^^^^^^^^^^^^^ 

characteristics  of  M30  propellant.  A  closed-form  expression  for  the  failure  time  (time  at  maxnnum 

stress  level)  is  obtained  by  setting  the  SPess-rate  in  (5)  equal  to  zero,  - 


dig,,):,::: 


di 


■+: 


2£ 


.  6  1/r. 
die  ) 

dt 


the  chain  rule  of  differentiation  for  the  first  term  in  (20)>  4g,^)/r/t  =  d(g,^)/JS,  d 
some  algebraic  manipulatibn,  the  princifiaf^m^^  maximum  stress,  % 

constant,  i.e.. 

S^iit,  and  after  ( 

at  failure  -  ^ 

where  3  1.5  a  r /(log, (1-n)).  Since  log, ^  (S^)/log„,  damage  parameter  at 

the  failure  stress  is  constant,  i.e.,  The  time  required  to  reach  maximum  stress,  or 

tciiltifc  tim&,  det^rttiitted  by  solving  (15)  for  subject  to  condition 

(22):;;:;; 

lit 

The  log  failure  stress  versus  log  failure  lime  is  plotted  in  Figure  6.  The  predictions  (solid  symbols) 
compare  well  with  observed  failure  stresses  and  observed  failure  times  (open  symbols).  The  failure^^ 
stress  is  also  observed  to  nionotonically  increase  with  strain  rate  (open  symbols  in  Figure  7)(  The 

increase  in  failure  stress  with  strain  rate  is  also  predicted  by  (23)  (solid  symbols)  and  the  I 

weak  strain  rate  dependence  is  given  by  o,  =  59  e  11,5  failure  stress  is  insensitive  to  strain 

rate  if  n  <  <  1  in  (22)  and  (23).  The  theory  predicts  that  the  strain  at  failure, 

-  tfC 

is  independent  of  strain  rate  in  materijUs  if  r  ^  nq  (Harper^^^^^ 

ihg  (22)  into  (24)  with  r  =  nqv  Tlhe  strain  rate  independence  of  failure  strain  is  alsa  predicte  the 
theory  in  niaterials  if  r  =  1  and  n<<  i;  as  in  linear  elastic  materials  with  dainagev)B 
conditions  are  approximately  true  for  M30  material  constants*  and  this  may  explain  why  measured 

yield  strains  (approximately  equal  to  failure  Strains)  in  M30  are  insensitive  to  strain  rate  (Gaz^^^ 

1991 ;  Gazonas  and  Ford  1992).  In  addition  to  being  strain  rate  independent,  the  failure  strain  iS  seen 
to  take  on  a  constant  value  since  the  slope  of  faimre  Stress  verSu  timei  curve  (Figure  6)  and 

the  Slope  failure  stress  versus  strain  rate  curve  (Figure  7)  are  equalin  magnitude)  but  opposite  in 

sign. 

stresses  are  determined  if  the  prodUctof  the  relaxation  modulus  and  the  strain  history  in  (6)  is  an 
integrable  function  Strain  histories  of  arbitrary  complexity  can  be  approximated  using  a  large 
number  of  constant  strain  rate  ramp  functions.  As  an  illustration^  an  input  history  that  consists  of 
two,  successive,  constant  strain  rate  ramps  is  expressed  by 

e(t)  =  Ej  t  H(t)  +  (E^  -  *ej )  t  H(t  - 1^  ,  (25) 

where  Ej  and  are  the  constant  strain  rates  of  each  ramp  and  t^  is  the  time  of  application  of  strain 
rate  U^  .  time-dependent  stresses  are  determined  by  substituting  (25)  into  (6).  Two  examples  that 
illustrate  such  ramp  inputs include:  1)  time-dependent  stress  predictions  for  Stress-relaxation  of  M30 
after  a  period  of  constant  strain  rate  straining  appears  in  Figure  8,  and  2)  a  “ballistic-like”  input, 
Simulated  by  a  concavemp,  two-ramp,  input  history  ,  with  the  strain  rate  in  the  second  ramp  an  order 
of  magnitude  greater  than  the  strain  rate  in  the  first  ramp  (Figure  9).  A  comparison  of  the  constitu¬ 
tive  response  of  M30  under  a  variety  of  input  histories,  such  as  those  given  in  this  section,  will 

provide  a  more  general  verification  of  the  theory.  Experimental  programs  should  include  input 
histories  that  are  Similar  to  those  experienced  by  the  propellant  in  the  gun  during  firing,  ^  ^ 


REMARKS 

A  nonlinear  theory  of  viscoelasticity  with  damage  has  been  shown  to  accurately  model  the 
constitiitive  response  of  M30  gun  propellant  in  ud  isothermal  compression.  The  exponential 
damage  function  originally  proposed  by  Schapery  (1981)  and  used  to  describe  uniaxial  deforma¬ 
tion  and  failure  in  ice  (Harper  1986)  was  found  unsuitable  for  describing  deformation  and  failure  in 

M30  propellant.  Instead,  ah  ''ellipsoidal'’  damage  function  was  determined,  directly  from  the  data;^ 

which  accurately  predicted  worksoftening  behavior  under  monotonically  increasing  deformation.  In 
addition,  time-dependent  predictions  of  stress  versus  time,  and  failure  stress  versus  failure  time,  and 
tailure  stress  versus  strain  rate*  quantitatively  agree  with  experimental  results  from  constant  strain 
rate  tests  on  the  propellant.  The  observed  insensitivity  of  failure  strains  to  strain  rate  in  M30 
(Gazottas  1991 ;  Gazohas  and  Ford  1992)  is  a  result  that  is  also  derivable  from  the  theory.  Future 
work  is  planned  to  verify  the  generality  of  the  constitutive  model  under  more  complex  input  histories 

such  as  those  described  in  the  previous  section.  However,  sighifiOaiif^^^^^^^^ 
observed  in  the  propellant  after  unloading  may  pose  difficulties  fofthe^m 
model  to  include  femperatute-dependent  behavior  with  reduced  time  vmiahles  is  also  currenUy  u 

development. 
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Table  1 


Axial  Sixers  versus  Tin^e  at  Various  Strain  Rates  in  M30  QUU  PfQpgU<ffl^ 


Component  % 

Nitrocellulose 

28.0 

%  NC  Nitration 

12.6 

Nitroglycerin 

22.0 

NitrOguanidine 

48.6 

Ethyl  Centralite 

2.0 

Length  (mm) 

25.4:;.:: ':■ 

Diameter  (mm) 

15 


Damage  (D  ) 


intentionally  left  blank. 
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9.  APPENDIX  A:  DERIVATION  OF  GENEMLIZE^  DAMAGE  FUNCTION 

::;i::^;::;::::AND:DAMAGRPAR 

The  derivation  of  the  gerieralized  damage  function  and  parameter  given  by  (  I 

(12)  in  the  tnain  lext  begins  by  taking  the  fet  Unte  defivadve  (dotted  quantity)  Of  the  function  f  in 
(3), 


and  collecting  ternivS  to  Obtain, 


Since  gj^  g^(a)  and  g,=g,(S^),  ordinary  derivatives  of  these  quantities  are* 


rfo  rfc 


”2  dS„  dt 


(A.3) 


(A,4) 


Substitution  of  (  A,3)  and  (A.4)  and  the  first  time  derivative  of  the  damage  parmeter  given  in  (4) 
from  the  main  text  into  (A.2)  yields  the  following  nonlinear  ordinary  differential  equation  in  a,  of 
Bernoulli  form, 


02 


where,  f,  =  fj(t),  A  change  of  variables  in  (A.5)  using. 


(A.5) 
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results  in  the  following  linear  ordinary  differential^e 


where 


u  +  p(t)u  s  h(t) 


£i 

f  r 


and  ii(t) 


r  g. 


An  integrating  factor  p(t)  for  (A.8)  is  obtained  by  letting, 


(A.6) 

(A.7) 


Iff 
r  J  ? 


5^  In  1  f  I  +  c 


(A.10) 


Substitution  of  (A.6)  into  (A:12),  using  the  identity  (10),  f  ^  e^y  lTom  the  niain  text,  and  letting 
=  1,  one  arrives  atthe  generalized  damage  function  (11)  and  damage  parameter  ^ 
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